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Abstract 

We show that the so called Robinson- Trautman-Maxwell equations 
do not constitute a well posed initial value problem. That is, the 
dependence of the solution on the initial data is not continuous in 
any norm built out from the initial data and a finite number of its 
derivatives. Thus, they can not be used to solve for solutions outside 
the analytic domain. 

PACS numbers 03.30.De, 04.40.Nk 

1 Introduction 



One of the most challenging problems in GR is to construct non station- 
ary space times for which a null boundary Scri can be attached. Those 
space times are particularly useful since the radiation fields at Scri repre- 
sent the gravitational and/or electromagnetic radiation emitted from some 
source[l][2][3][4]. Finding explicit solutions of the Einstein equations with 
this property, however, is a formidable task unless some extra conditions are 
assumed for the space time under consideration. Sometimes the set of con- 
ditions simplifies the problem in such a way that either and explicit solution 
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can be found or the resulting equations are sufficiently simple that they can 
be integrated numerically. 

In particular, if we assume that the space time admits a shear free and 
twist free null congruence, the field equations adopt an extremely simple 
form. In 1960, Robinson and Trautman (RT from now on) obtained a class 
of Ricci flat space times where all the field equations but one could be ex- 
plicitly integrated [5]. This last equation gives the dynamical evolution of the 
radiating gravitational field at Scri and the solution is regarded as a radiating 
space time since it has a non trivial Bondi mass that obeys the Bondi mass 
loss equation and decays exponentially to a stationary space time. Although 
one could argue that gravitational waves should not decay exponentially but 
rather with a power law formula, the RT solutions have been used as a work- 
ing model and several stability theorems have been proved to show that it 
might be regarded as a physical solution[6][7][8][9]. 

It is an interesting problem to generalize the RT equations to include 
electromagnetic radiation and call Robinson- Trautman-Maxwell (RTM from 
now on) space times to the new class of space times. In this note we address 
this problem, we write the RTM equations and then show that if the total 
charge is different from zero the field equations are unstable against linear 
perturbations, but worst than that, the growth rate increases without bound 
as the wave number of the perturbation tends to infinity. Thus, by a simple 
argument we conclude that the RTM equations are not well posed and so 
they can not be used as an ansatz to construct physical solutions. 

In Section 2 we write down the relevant field equations whereas in Sec- 
tion 3 we analyze the stability of those equations. The use of the stability 
theorems and its implications on the existence of those space times is given 
in the Conclusions. 

2 The equations 

The type II algebraically special space times have two principal null direc- 
tions. Thus, it is possible to find a null tetrad such that two of the scalars 
built out of the Weyl tensor by contracting it with the frame, usually denoted 
by ip , and ipi, vanish. It can be shown that all the radial equation of the 
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Einstein Maxwell field equations can be integrated leaving just the angular 
and time evolution to be solved. 

If in addition we impose the extra condition that the tetrad is shear and 
twist free the resulting equations adopt a very simple form when written 
in a suitable coordinate system instead of the standard Bondi coordinates 
(u, q, <T)[10]. This particular coordinate system (r, q, q) is called Newman 
Unti and the relationship between the NU and Bondi time r = T(u,q,q) is 
one of the basic variables for the problem. It can be shown that the Ein- 
stein Maxwell field equations constitute a set of 4 equations for the variables 
(0i, <p2, V,ip2) where <pi and 4>2 represent the non vanishing Maxwell scalars, 
^ 2 is the Weyl scalar thar represents the mass aspect, and V = X'with 
u = X(r,q,q) the inverse function of T. With this in mind we write the 
RTM eqs. and suitable definitions below. 

Maxwell 1. Defining q := (f)\V 2 we write the first Maxwell equation as 

9g = =>- q = q(r) 
Restricting the freedom left in r we can set q = const. 

Maxwell 2. We define $ := 4> 2 V and the second Maxwell equation reads 

5$ = -[qV- 2 ]', 

from which we obtain an evolution equation for V, 

2qV' = V 3 m. 

GR 1. Starting with 

v- 3 d(^ 2 v 3 ) = 2x0!0 2 , 

using 0! = prjand defining \ := — faV 3 , we can rewrite the above equa- 
tion as 

dx = -2xg$ 

An equivalent eq. is given by 

4xq 2 V + V 3 (V 2 x) = 

with V 2 x = S9x the Laplacian operator on the unit sphere. 



3 



GR 2. 

(v~ 3 xY - {(39) V + 233V - v _1 9V • gV} + xv> 2 2 = o 

which can be rewritten as 

4xgY + 3x^ 2 (V 2 x) - 4xg 2 (V 3 (VV + 2VV) - V 2 (VV) 2 ) + ( V x ) 2 = 
Thus the non trivial eqs. to solve are 

4xgV + V 3 (V 2 x) = (1) 
4xgY + 3x^ 2 (V 2 x) - 4xg 2 (V 3 (VV + 2VV) - V 2 (VV) 2 ) + ( Vx) 2 = 

3 The stability of the equations 

To show that the equations are not well posed we study their linearized 
version, first off a constant solution, which allows a full treatment and then 
around any solution, using a theorem by Strang [[11]]. 

3.1 Lack of well posedness around constant solutions 

To analyse the stability of the above set we first write down their linearized 
versions off the Reissner- Nordstrom solution (given by q = const., X = Xo = 
const., V = 1), 

4xg 2 x' + 3xo(V 2 x) - 4xg 2 (VV + 2VV) = 

4xgV + (V 2 x) = 
Taking a time derivative on the top eq. gives 

4xg 2 x" + 3xo (VY) + (V 6 x + 2V 4 x) = 
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Looking for a solution proportional to an eigenfunction of V 2 with eigen- 
value k 2 we get 

4xq 2 X "-3Xok 2 x'-(k (i -2k 4 )x = 
and so the solutions have a time dependence of the form e at , 
Axq 2 a 2 - 3xok 2 a - (k 6 - 2k A ) = 0. 

with a given by, 

_ 3 X ok 2 ± i/{3xok 2 ) 2 + 16 xq 2 {k 6 -2k 4 ) 

For fixed values of k and arbitrary values of t, either positive or negative, 
one branch always blows up and the linear system obtained studying per- 
turbations around a constant solution is unstable. Worst than that, for any 
fixed value of t , as k — * oo , both roots of a blow up, i.e., the growth rate 
of the perturbation increases without bound as the wave number increases. 
This implies that the system is not well posed in the sense that the solution 
is not a continuous function of the initial data when the topology is given by 
a norm controlling a finite number of derivatives. 

To see this, take a sequence of initial data of the form 

Xk = j^fk(x),x' k = °> 

where fk(x) an eigenfunction of V 2 with eigenvalue k 2 , and p is the maximum 
number of derivatives controlled by the norm. Thus, as k — > oo the norm 
of the initial data remains bounded, on the other hand, since the solution 
corresponding to this data is of the form fk{x)(e a+t — e a ~') the norm at time 
T > would have grown by a factor e a+T (or by a factor e OL " T if T < 0). 
In both cases the norm of the solution at T would grow with no bound as 
k — > oo showing that there can not be a bound of the solution in terms of 
the bound on the initial data for any finite time and for all initial data. 

This result is not only valid for linear perturbations but, as we shall see 
below, extends around arbitrary smooth solutions. The reason is that by 
being essentially a high wave number phenomena we can consider an arbi- 
trary small neighborhood of an arbitrary point in space time and localize 
the perturbation there. In such a small neighborhood the background solu- 
tion can be considered of constant coefficients and so the previous analysis is 
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valid. The perturbation need not be of the order of the background solution 
to show ill posedness, for it is enough to see that the difference of solutions 
growth as in the linear case and that difference can be made arbitrarily small 
on the initial data. 

Even if we find a real analytic solution of the equations any C°° arbitrarily 
small infinitely different iable perturbations of its initial data will produce 
arbitrary large solutions when we evolve the data with the field equations. 

We turn now to the more general analysis of the stability around arbitrary 
smooth solutions. 

3.2 Lack of well posedness around an arbitrary solu- 
tion 

To study the general case we invoque a theorem of Strang [11] asserting that 
if a linear, smooth coefficient system is well posed in the L 2 norm of its 
components, then the principal part of it with its coefficients frozen at any 
point must also be a (constant coefficient) well posed system. This theorem 
illustrates the fact that the issue of well posedness is a microlocal problem 
or, equivalently, a high frequency phenomena. We shall use Strang's theorem 
in the double false way, that is, showing that if the principal part is not well 
posed, then neither is the variable coefficient system. 

In order to apply Strang's theorem we first take the gradient of the first 
equation of the full system (1), and then consider the principal part of that 
new system, which, for quasilinear systems, is the same as the one corre- 
sponding to their linearization, at an arbitrary point p: 

x' - V 3 (p)V 2 V ■ V = 0, 

where we have defined the new variable V = VV, and so the vector of 
variables is now (V, V, x)- 

The eigenvalues of this system are 

e 

a ~ 2 v ^qV 3 (pY 

that is, the limit for large k of the ones previously found for linearizations 
off constant solutions with V — 1. We clearly see that the frozen coefficient 
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system is ill posed and so, by Strang's theorem, is the linearized equation of 
the full system around any smooth solution. 

3.3 The RT branch 

It is interesting to see why the RT equations can be recovered in the limit 
when the charge q tends to zero, and so why are they well posed in the 
negative time direction. To see that we now assume that the solution can be 
written as 

x = xo + q 2 x 

Thus the linearized eqs. for \ and V read 

4xg 2 x' + 3xo(V 2 x) - 4x(VV + 2VV) = 

AhV + (V 2 x) = 
Taking a time derivative on the top eq. gives 

4xg 2 x" + 3 X o(V 2 xO + (V 6 * + 2V 4 x) = 

The Fourier transformed eq. is 

4Wx"-3Xok 2 x'-(k e -2k 4 )x = 
the Laplace transform gives 

4xgV - 3xok 2 a - (k 6 - 2k 4 ) = 

so 

_ 3 Xo k 2 ± f/{3 X ok 2 ) 2 + 16xq 2 (k e - 2k 4 ) 
a ~ 8^ ' 

i.e., we obtain the same roots as before. Taking xo < 0, and the positive 
root and then the limit q — * yields the standard RT linearized solutions. 
Thus, in the limit, we see that there is a branch whose growth rate is bounded, 
in this case for positive times. Those are the ones which in that limit go to 
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the usual RT linearized solutions. Even the structure of the RT solutions 
is puzzling, for one can only evolve them to the future but not to the past, 
while Einstein's equations in their symmetric hyperbolic formulations allow 
for evolution both into the past and future. 1 



4 Conclusions 

We have shown that the RTM equations are not well posed and so they can 
not be used as an ansatz to construct physical solutions beyond the analytic 
domain. 

This is rather unexpected, since it is well known that the full set of 
Einstein's equations without any special conditions is a well posed problem 
under suitable gauge conditions. The problem seems to arise from the special 
condition assumed for the space time which forces some Weyl scalars to 
vanish throughout the evolution. This condition is not a typical local gauge 
condition on the initial data set. Thus, even if we can find a real analytic data 
for which a solution exists, an arbitrary C°° perturbation would not preserve 
the algebraic condition and that would be reflected on the appearance of a 
unbound growth. 

This lack of well posedness makes the system totally unsuited for finding 
numerical solutions. Any numerical scheme is prone to errors (particularly 
truncation errors) and even if we start with analytic data, the errors (consid- 
ered as perturbations) would completely overcame the exact solution. There 
won't be any way of having a stable converging algorithm. 

Moreover, the above analysis shows that if we perturbe a Robinson Traut- 
man space time adding a small amount of electric charge, the whole construc- 
tion blows up. Since we cannot prevent the presence of tiny amounts of excess 
charge in any compact source it appears that RT space times do not represent 
physical sources. 



^n some cases, like the one where the initial surface extends to future null infinity, 
only future evolution is possible for the maximal domain of dependence, but even in this 
case local solutions on the causal domain of the initial data surface is possible. 
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